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Recent laboratory data (Durham et al,, 1983; Kirby et al, 1985) reveals that H 2 O ice 
has a strongly non- Newtonian (power-law) rheology at temperatures appropriate to 
the near-surface regions of the Galilean satellites i.e. that 

i = Aexv(-Q/RT)T n (l) 

where e is strain-rate, r is effective stress and A — 2.51 X 10~ 26 Pa“ 3,7 s -1 , Q = 
27 kJ mole -1 and n = 3.7 for 158 K < T < 195 K. Numerical calculations in- 
corporating this rheological law indicate that crater relaxation in a non-Newtonian 
rheology has significantly different characteristics to viscous relaxation in a Newtonian 
medium (Thomas and Schubert, 1986). These arise principally by virtue of the fact that 
viscosity rj varies as r\ ~ r 1-n . Thus high stresses are associated with low viscosities, 
which permit relatively rapid relaxation flow. However, once the initial stress is re- 
duced, viscosity increases and relaxation flow becomes markedly slower. Consequences 
of this phenomenon observed for numerical simulations of crater relaxation include the 
formation of a high-stress, low- viscosity region immediately below the center of the 
crater; a relatively rapid relaxation rate intially, followed by a marked slowing of the 
relaxation flow as stresses decrease (and viscosities increase) and, concomitantly, much 
greater differences between relaxation timescales for small vs. large craters than would 
be expected for a Newtonian medium. 

Calculations carried out for craters of varying sizes in a medium with constant 
temperature T = 173 K yield values for the crater relaxation time t e (defined as the 
time required for the crater depth to become 1/e of its original value) that appear to be 
too small to account for the observed retention of craters on Ganymede and the other 
icy satellites. Typical values were t e fn 10 2 y for craters with a diameter of D = 300 km 
and t e « 10 5 y for D — 50 km. Although the rate of crater relaxation at times greater 
than t e is much slower than the initial rate, as mentioned above, as many as IQ 6 en- 
folding times are available to reduce the remaining topography. Such a calculation is 
thus seriously in conflict with the observed crater population of the surfaces of the icy 
satellites. 

In an attempt to reconcile this conflict, we consider possible explanations for 
the much slower relaxation rate of craters on the icy satellites. It is possible that 
an admixture of silicates in the surface ice regions of the icy satellites may raise the 
viscosity to some extent. However, even a silicate volume fraction of the order of 
0.6, implausible for the icy satellites, whose surfaces seem to be predominantely ice, 
only increases the viscosity (and thus the relaxation time) by at most two orders of 
magnitude (Friedson and Stevenson, 1983). 

A rheological transition to a Newtonian law (where n = l) will not occur 
unless the the Newtonian rheology permits a greater strain-rate for a given stress 
than the non-Newtonian rheology. One cannot thus increase viscosity by appealing 
to a rheological transition mechanism. In fact, the critical stress r c below which a 
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Newtonian rheology is dominant, given currently accepted values for Newtonian ice 
viscosity laws (Q = 57.7 kJ mole -1 ) and ( A = 9.06 x 10 -16 Pa“ 1 s -1 ) (Reynolds and 
Cassen, 1979) is of the order of 10 Pa, many orders of magnitude less than the stresses 
associated with even small craters. 

It may be argued that the grain size of the ice in the upper layers of the 
icy satellites is unknown, and this may render inapplicable the experimental results of 
Kirby et al. (1986). While the creep of polycrytalline ice is sensitive to the grain size 
at temperatures close to the melting point (Baker, 1978; Goodman et al., 1981), it is 
probably of less importance at low temperatures (S.H. Kirby, personal communication, 
1986). In any event, this effect alone is not likely to be responsible for the orders of 
magnitude increase required in viscosity. 

A solution to the problem can be achieved by accounting for temperature 
variations in the surface regions of the icy satellites, and utilizing more recent exper- 
imental data (Kirby et al., 1986). This data, the product of an over 200% increase 
in experimental runs at low temperature, yields n = 4.8, A = 2.10 x 10 -33 Pa _4-8 s -1 
and Q = 29 kJ mole -1 . The significantly higher value of n is observed to enhance the 
non-Newtonian phenomena associated with crater relaxation and mentioned above. 

To incorporate the effects of a temperature gradient, we take the case of 
Ganymede (T 3Ur f ace = 120 K) and assume a temperature gradient dT jdz = 1 K km -1 , 
consistent, for the present day Ganymede, with global thermal models (Schubert et 
al., 1981). These models predict the onset of convective motion at a depth where 
T ss 0.6 Tm = 164 K. The adiabatic temperature gradient for the icy satellites is small, 
and thus we assume that the convecting region is isothermal, with T = T m . Thus 
the maximum temperature determining the viscosity for crater relaxation is 164 K, 
although higher temperatures may be appropriate to the early history of Ganymede. 
We further assume that the dominant viscosity is that determined by the temperature 
T a at a depth z = D/2 (which lies in the center of the relaxation flow region (Thomas, 
1986)) and a typical stress r a determined from numerical analysis of craterform stresses 
(Thomas and Schubert, 1986). We can thus construct estimates of relaxation times for 
various craters as follows: 
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Kirby et ah (1986) did not observe viscous behavior for T < 158 K, and thus viscous 
relaxation is not expected for the craters of smallest sizes. Relaxation times are cal- 
culated for the D = 75 km crater, a borderline case for which viscous relaxation may 
occur. 

These results are in good qualitative agreement with observations of the sur- 
face of Ganymede: they account for the paucity of well-preserved large ( D > 100 km) 
craters without requiring any ad hoc mechanisms such as impact melt-lubrication 
(Croft, 1983) or acoustic fluidization (Melosh, 1979) while allowing for the long-term 
retention of smaller craters. Further research will concentrate on the quantitative 
comparison of calculated relaxation times to the observed population of craters on 
Ganymede. 
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